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Abstract 

In a fully relativistic approach, a RLSM description of nuclei colliding at ultra-relativistic energies can be 
formulated within the framework of a classical transport theory. The valence quarks of the nucleons are 
described through collections of classical point-like particles moving in the continuum. They are coupled to 
soft gluon fields which are described through the Yang Mills equations on a gauge lattice. In a first step, 
we focus on the range of low-pt interactions. Results from numerical model simulations of pure gluonic 
nucleus-nucleus collisions on SU(2) gauge lattices in 3+1 dimensions are presented. They show an effect 
which we call the glue burst. 



The state of super-dense nuclear matter called the quark gluon plasma is expected to occur in the central 
kinematic region of ultra-relativistic heavy-ion collisions. In view of upcoming experiments at the RHIC 
in Brookhaven and the LHC at CERN, it is one of the most challenging topics in the physics of ultra- 
relativistic heavy-ion collisions to develop a coherent description of the formation of the quark gluon plasma 
state. Many of the descriptions which have been developed over the past years are based on the idea of 
a perturbative scattering of partons within transport models. One of the problems in these descriptions 
concerns the initial state of the colliding nuclei. The transport equations start from probability distributions 
of partons in the phase space. In reality however, the states of the nuclei are described by coherent parton 
wave functions. The incoherent parton description becomes inadequate at exchanges of small transverse 
momenta. A four years ago McLerran and Venugopalan proposed [1] that the proper solution of these 
difficulties is the perturbative expansion not around the empty QCD vacuum but around a vacuum of the 
mean color fields which accompany the quarks in the colliding nuclei. This idea motivates a combination 
of the parton cascade model [2] with a mean field description of the color fields. Before this can be done 
successfully, it is important to study the non-perturbative dynamics of the color mean-fields themselves. 
In Fig. 1, a typical scenario of a collision is displayed for times shortly before and after the collision. The 
figure depicts the basic idea behind our model of a combined particle-gauge lattice system. The collision 
describes a scenario in which two media are involved: The valence quarks, described through collections of 
point-like color charged particles (solid lines) and the gluon fields (helical lines). While the quarks interact 
essentially only through the exchange of gluons, the gluon fields are self-coupled. 
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The right-hand part of Fig. 1 displays the basic concept of our gauge lattice transport model [3] . The soft 
color fields are represented on a SU(2) gauge lattice in three dimensional Euclidean space and are described 
by the Yang-Mills equations. The quarks are represented as collections of color charged particles moving in 
the continuous background in position space. Their trajectories are solutions of the Wong equations. The 
particles are coupled with the lattice by energy and color exchange mechanisms. In SU(2) a nucleon initially 
is represented through a pair of quarks with opposite momenta and opposite color charge. In principle the 
model allows the inclusion of collision terms to describe the hard collisions between the particles. This idea 
has recently been submitted [3] . 

In recent simulations of a central collision of two nuclei with the above described model [4] it turned out 
that in addition to large transverse gluon radiation there exists a contribution to the transverse dynamics 
from "soft glue field scattering". This leads to the question, how important is really the mean field dynamics 
for the transverse dynamics. Due to non-linear effects, the gluon mean field might contribute more to the 
transverse momenta than naively expected. We, therefore, go first one step back and focus on the pure gluon 
mean field dynamics in the non-perturbative regime and leave out the particles. 

In short denotation the homogeneous Yang-Mills equations read 2? M JT Ml/ = 0. An lattice version of the 
continuum Yang Mills equations is constructed by expressing the field amplitudes as elements of the corre- 
sponding Lie algebra, i.e. £s,k,Bx,k G LSU(2) at each lattice site x. We choose the temporal gauge Ao = 
and define the following variables. 

Ug t i = exp(— igaiAi{x)) = Ut +l _ t Ug.ki =U3,kMx+k,iMx+k+i,-kMx+i,-i (1) 

In adjoint representation the color electric and color magnetic fields are expressed in terms of the above 
defined link variables Ux,i and plaquette variables Ux,ki in the following way. 

= 7~T^ ' Bs,j = , tjki (Us,m — (2) 

igiij ^yukiii 

The lattice constant in the spatial direction I is denoted by aj. We choose Us,i and Ex,i as the basic dynamic 
field variables and numerically solve the following equations of motion. 

iis,k{t) = iga k Sg }k (t)Ug^(t) (3) 

3 

isAt) = 2ga ia2 a z ~ U Im® ~ U l-i^ fa-iM® - U%_ lM (t)) M g _,,,(0 }■ (4) 

We first study the collision of plane wave packets which have initially constant amplitude in the planes 
transverse to the wave vector. We choose periodic boundary conditions implying the topology of a 3-torus 
for our lattice. The size is chosen as 10x10x2000 lattice points. For a collision, we initially arrange two 
Gaussian wave packets with average momenta (0, 0, ±fc z ) and width Afc; between the maximal and minimal 
momenta on the lattice: k™ m <C Afc; <JC k t <Si k™ ax The polarization in color space is defined through the 
unit vectors n c L for the left and n c R for the left (L) and right (R) moving wave packet, respectively. 
We express the initial conditions through the gauge fields [5] 

Ar /l = S n n c L/R 4>(=ft,z-Z R/L ) (5) 

where the scalar factor <j>(t, z) defines the initial wave packet 

<f>(t,x 3 ) := <j) exp(-]-Akl(t + x :i ) 2 )cos(k3(t + X3)). (6) 



The amplitude factor is defined as cf) := \J (2Afc 3 ) / (y^Kaks) where the parameter a denotes the transverse 
area per quantum contained in the wave field. We choose the polarizations in color space Hr — ul = (0, 0, 1) 
for the Abelian case and — (0,0, 1), ul = (0, 1,0) for the non-Abelian case. Once the initial fields are 
mapped on the lattice, the time evolution of the collision starts from the superposed initial conditions 

7/(0) _ 7/(0) 7/(0) c(0) _ c(0) , c(0) (7 s 

M S,i — M 3,i,L ' M SJ,.R> °x,l — °x,LR "•" °3,l,L> (') 

respectively. A linear superposition of solutions obeys the Yang-Mills equations only if these solutions have 
overlap zero. Therefore, the initial separation AZ of the wave packets may not be chosen to small. 
The calculation contains four parameters. The relative isospin polarization which is parameterized through 
the angle 9c- The average momentum of the wave packets k and their width A&3 and the coupling constant 
g which can be rewritten in terms of the parameter a as g' = g/y/a by simultaneously rescaling the field A 
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as A' — \J~cfA. Consequently, the system shows the same dynamics for different values of g and a, as long 
as the ratio g' = gjyfa is kept fixed. 

In the following we present results from a simulated collision in the non-Abelian case with the parameters 
ks — 7r/(2o3), A&3 = 7r/(100a3), g = 1, a = 1. Further, we used a; = 0.1 fm and At = 0.01 fm. Subsequently, 
we refer to the direction of the collision axis (the z-axis) as the "longitudinal direction" and to directions 
perpendicular to the collision axis as the "transverse directions". Accordingly, we define the transverse 
and longitudinal energy densities of the color electric field and the transverse and longitudinal energy flow 
distributions 
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w^\t,z) = y^aia2 /^Trfg^jgg,;), w { jf\t, z) — aia2Tr(6 Si3 £ Si3 j (8) 

x,y i — 1 x,y 

Fig. 2 displays {fin, z) plotted over the z-coordinate at various time steps t n as indicated on top of the 
curves. 
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Figure 1: Transverse energy densities WT{t,x^) are displayed for selected time steps t n 

At initial time to the curves {to, z) of the wave packets are completely filled indicating the strong 
oscillations according to k$ — 77/(203). After 1700 time steps the wave packets are colliding and have reached 
maximum overlap. The smooth white zone at the bottom indicates the appearance of long wavelength modes 
in the Fourier spectrum of the common field. These low frequent modes have been excited in the penetration 
time. At time step £6000, almost all the energy which has originally been carried by short wavelength modes 
around A?3 is transmitted into long wavelength modes. This can be seen from the high frequent oscillations 
on the surface of the big white zones in the two receding humps. At time step isooo however it seems that 
energy is transmitted back into high frequent modes and finally the wave packets start to decay around time 
step tioooo- Simultaneously, we display the corresponding longitudinal energy densities w^\t, z) in Fig. 3 
for the same time steps t n but on a logarithmic scale. 

We remember that the wave packets were initially polarized into the transverse x-direction and consequently 
has to be zero as long as the wave packets propagate freely. However, when the two colliding wave 
packets of different color start overlapping, we observe an increasing longitudinal energy density in the 
overlap region around the center of collision at z = 100 fm. Fig. 3 clearly shows that w ® {t, z) grows rather 
fast from time step iiooo until time step ii7oo where it reaches a maximum and has grown exponentially by 
more than two orders of magnitude. For larger times the hump decreases and practically disappears at £3000. 
Around tgooo, however, the longitudinal energy density grows again at the positions of the receding wave 
packets. After 10000 time steps w^ 1 (z) has increased by five orders of magnitude. Now the question arises 
whether the energy deposit in longitudinal links can be associated with fields propagating into transverse 
directions. The transverse and longitudinal energy flow densities displayed in Fig. 4 and Fig. 5 are defined 
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A direct comparison of Fig. 4 with Fig. 3 shows that the transverse energy flow distribution st(£, z) behaves 
very similar as compared to the longitudinal energy distributions. We call the phenomenon resulting in the 
strong increase of the transverse energy flow at large times the "glue burst". The most interesting feature is 
the time delay extending from time step ti600 up to teooo- Repeating the calculation for increasing g- values, 
we found that it scales with g~ 2 . The explanation can be found through a carefull analysis of the Yang-Mills 
equations. The growing peak in the center of collision clearly shows that the wave packets interact strongly 
due to the non-linearity of the Yang Mills equations. The wave packets, however, are restored after they 
have passed through each other and look for a while almost unchanged. Suddenly, around time step teooo 
they start to decay. At time step tioooo the transverse energy flow density has grown by three orders of 
magnitude in comparison to the one at time step teooo ■ Afterwards the energy flow density shows no further 
increase. A calculation for the Abelian case with same parameters shows no fragmentation of the wave 
packets, no transverse energy currents and no peak-like structure in the overlap region. 
Fig. 5 displays the corresponding longitudinal energy current densities SL(t,z). As a function of time they 
behave similarly as the transverse energy densities wjf '(t, z) shown in Fig. 2. The graphs in Fig. 4 do 
not appear to be completely symmetrical. Matinyan et al. have shown, that the Yang Mills equations are 
sensitive to small perturbations. The difference in shape of the right moving and left moving energy flow 
distribution reveals a small difference in the initial conditions by reading them slightly asymmetrically on 
the lattice. This effect is enhanced in our example because we have chosen a very short initial wavelength 
&3 = 7r/(2o3). At high oscillations in position space, slight phase shifts can strongly change the shape of the 
density distribution. Nevertheless, the total integrated energies of both wave packets agree exactly. Fig. 2 
to Fig. 5 are created from data obtained in one run of our computer code. A comparison of Fig. 5 with 
Fig. 2 shows that the asymmetry comes from the B-fields which seem to be more sensitive to perturbations. 
Further, a comparison of Fig. 4 with Fig. 5 shows that SL{t,z) in Fig. 5 decreases when ST(t,xz) in Fig. 4 
increases. At time step ii4ooo a sizable fraction of the longitudinal energy flow is converted into transverse 
energy flow. We have repeated the calculations at k^ — 7r/(100a3) and found qualitatively the same results. 
Here, however we demonstrate that our results are even found at the critical average longitudinal wavelength 
of 4a3 which is just twice as long as the shortest wavelength of 2a2 at the lattice cut-off. We also observe 
that the transverse energy current increases for larger k$. 

In a study of colliding wave packets of finite extension into one transverse direction we use the same initial 
conditions as in the case of plane waves but multiply the function <j> in Eq. (6) with a factor exp (— (Ak y y) x ). 
In the simulation of the collision shown in Fig. 6 below we used the parameter values Ak y = 1/(40(22) and 
A = 8 on a 4x200x800 lattice. The coupling constant was taken g = 6.0 in order to reduce the time delay. 
The large value of A leads to a width of 8 fm in the y-direction leaving a space of 6 fm on each side of the 
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initial wave packet. Similar results as for colliding plane waves were found. The induced transverse energy 
currents, however, were much larger. Fig. 6 shows the transverse energy distribution at the time steps £100, 
£1500, £3500, £4000, £5000, and £6000. The glue burst occurs around time step £4000. The full explanation 
of this phenomenon will be published elsewhere. 

Note: The Fig. 6 exceeds the size accepted at the lanl-server and was therefore omitted. We recommend to 
download the postscript file from the web page http://www.phy.duke.edu/" poeschl/ under "Rsearch Related 
Links" and "Preprints". 
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